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Abstract 

We study gluon radiation in top quark production above threshold at high energy 
e + e~ colliders. We allow for the top quarks to be off-shell, considering radiation in 
both the top production and decay processes simultaneously. Our calculation includes 
all top width effects, spin correlations, and b quark mass effects. We study the effects 
of gluon radiation on top mass reconstuction and examine the interference between 
production- and decay-stage radiation, which can be sensitive to the value of the top 
quark decay width. 
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1 Introduction 



By virtue of its large mass, the top quark is unique. Because it is so heavy, the top quark 
decays before it can form hadrons |]J and its spin information is passed along to its decay 
products. Top's mass is comparable to the electroweak symmetry breaking scale, and its 
Yukawa coupling is suspiciously close to unity. All of this means that top physics will be 
an important and interesting subject in the near future, and the top quark may lead us to 
physics beyond the Standard Model. 

Future high energy lepton colliders — e + e~ and - can provide relatively clean 

environments in which to study top quark physics. Although top production cross sections 
are likely to be lower at these machines than at hadron colliders, the color-singlet initial 
states give lepton machines some advantages. Furthermore, the fact that the laboratory 
and hard process center-of-mass frames coincide greatly simplifies the reconstruction of final 
states. In addition, many of the top quark's couplings, especially those to the photon and 
Z° boson, can be easily studied there. 

The potential for precision studies of top physics at such colliders requires precision 
predictions from the theory, beyond leading order in perturbation theory. In particular, QCD 
corrections must be taken into account. Jets from radiated gluons can be indistinguishable 
from quark jets, complicating identification of top quark events from reconstruction of top's 
decay products. To make matters worse, emission may occur in either the top production 
or decay processes, so that radiated gluons may or may not themselves be products of the 
decay. Subsequent mass measurements can be degraded, not only from misidentification of 
jets but also from subtle effects such as jet broadening when gluons are emitted near other 
partons. 

In this paper we study the effects of real gluons radiated in top quark production and 
decay at e + e~ colliders. We consider collision energies well above the top pair production 
threshold, so although for definiteness we will refer to electrons in the initial state, our 
parton-level results apply equally well to collisions at the same energy. We allow for 

the top quarks to be off-shell, keeping the full width-dependent top propagator and retaining 
all spin correlations. Past treatments of QCD corrections have treated the top production 
@ or decay || processes separately, or in one case f| combined the two processes, but 
kept the on-shell (narrow-width) approximation. Gluon radiation for off-shell top has been 
treated previously in the soft gluon approximation || §]. Here we give an exact treatment for 
arbitrary gluon energies. We study properties of the emitted gluons, top mass reconstruction, 
and effects of interference between production- and decay-stage gluons that can be sensitive 
to the top quark width. Computation of the virtual QCD corrections for off-shell top is 
currently in progress JF| and will be combined with the present results for a full QCD NLO 
treatment of top quark production and decay. 

This paper is organized as follows. In Section 2 we discuss the matrix element calculation 
and Monte Carlo integration. We present numerical results in Section 3, focusing on gluon 
distributions, mass reconstruction, and interference effects. We present our conclusions in 
Section 4. An Appendix includes some comments about gauge invariance. 
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Figure 1: Feynman diagrams for gluon emission in top production and decay at lepton 
colliders. 



2 Cross Section Calculation 

2.1 Computation of the matrix element 

We compute the cross section for real gluon emission in top quark production and decay: 

e + e - _> 7 * ; z* -> ti(g) -> 6Vl/ + 6iy-^ . (1) 

At e + e~ colliders, no gluons are radiated from the color-singlet initial state. Final-state 
gluon emission can occur in both the production and decay processes, with gluons emitted 
from the top or bottom quarks (or antiquarks), as shown in Figure [l]. Emission from the 
top quark contributes to both production- and decay-stage radiation, depending on when 
the top quark goes on shell. Emission from the b quarks contributes to decay-stage radiation 
only. The separation of these contributions will be discussed below. 

We compute the exact matrix elements for the diagrams shown in Fig. [I] with all spin 
correlations and the bottom mass included, using the helicity methods of Kleiss and Stirling 
||. Working at the matrix element, rather than the matrix element squared, level has 
the usual advantages of numerical efficiency, and in our case has the additional advantage 
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that we can identify individual contributions and interferences between them. The explicit 
expressions for the matrix elements are complicated and not particularly illuminating, so we 
do not reproduce them heref]. 

We do not assume the top quark to be on shell; therefore we keep the finite top width T t 
in the top quark propagator and include all interferences between diagrams. We use exact 
kinematics in all parts of the calculation. We do not include radiation from hadronic W 
decays; the W bosons are assumed to decay leptonically and we integrate over their decay 
products in the results presented here. In practice radiative hadronic W decays could be 
largely eliminated through invariant mass cuts, but a detailed study is required to determine 
exactly how effective this would be. 

Some comments about gauge invariance are in order here. In the on-shell, or narrow- 
width, approximation, the top quarks are always on mass shell and the width is identically 
zero. In that case gauge invariance requires only diagrams with two intermediate top quarks. 
Once the top quarks are allowed to have a finite width and be off-shell, however, nonresonant 
diagrams — those with the same final state but only one or no intermediate top quarks — can 
also contribute. In the region of interest for top physics, viz., the region where the top quarks 
are nearly on shell, the poles in the doubly-resonant diagrams (those with two tops) cause 
their contributions to dominate. This is the motivation for the "double-pole approximation" 
in electroweak radiative corrections to W pair production |J. In the present case we use 
invariant mass cuts (see below) to guarantee that we are in a region of the phase space where 
the doubly resonant diagrams dominate the cross section. We can restore gauge invariance 
using the approximation described in the appendix. This gauge invariant result differs from 
the exact one by non-doubly-resonant terms, and this difference is very small numerically in 
the regions of interest. 

2.2 Production— decay decomposition 

As mentioned above, calculating at the amplitude level allows us to identify contributions 
from individual processes and their interferences. We are particularly interested in distin- 
guishing between contributions from gluons radiated in the top quark production and decay 
stages. This is directly related to reconstruction of the top quark momentum from its decay 
products, which in an experiment allows us both to identify top events and to measure mt- 
The presence of gluon radiation complicates the reconstruction because the emitted gluon 
may or may not be part of the top decay. If the gluon is not part of the decay, then it is 
represents a correction to top production and should not be included in the top momentum 
reconstruction: 

m 2 t « v\ = (Pb + Pwf = vlw ■ (2) 

If on the other hand the gluon is part of the decay, then it should be included in top 
reconstruction: 

™ 2 t ~ Pt = {Pb + PW+ Pgf = Plwg ■ ( 3 ) 

Being able to make this distinction turns out to be useful for purposes of efficient phase-space 
integration as well. 

3 A FORTRAN program containing the matrix elements can be obtained from the authors. 
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Although this production-decay distinction cannot be made absolutely in an experiment^ 
the various contributions can be separated in the calculation. For radiation from the b and 
b quarks, the assignment is easy: these contributions, corresponding to the two right-hand 
diagrams in Fig. |l], are clearly part of the top quark decay. However, as noted above, gluon 
emission from the top quark (or antiquark) contributes to both the production and decay 
stages; which is which depends on whether the top was closer to its mass shell before or after 
emitting the gluon. This condition corresponds to which of the two propagators from the 
top that emitted the gluon is numerically larger. 

We can make the separation in our calculation as follows . For definiteness, we consider 
gluon emission from the top quark, shown in the upper left diagram in Fig. [TJ. The matrix 
element for this diagram contains propagators for the top quark both before and after it 
radiates the gluon. The matrix element therefore contains the factors 

M « (-2 \ —tA \~2 1 —tA • ( 4 ) 

\Pwbg ~ m t + im t T t ) \p Wh - mf + vm t Y t ) 

The right-hand side can be rearranged to give 

M OC \ ; 5 \ ; — ] . (5) 

2>PWb * PWbg \Pwb ~ m t + lm t^t Pwbg ~ m t + lm J-t ) 

This separates the production and decay contributions to the matrix element. The first term 
in parentheses contains a propagator that peaks when p\ Vb = ntf, which corresponds to the 
condition for production stage, as in Eq. |2|. Similarly, the second term peaks for v\vbg = m t > 
which corresponds to decay emission as in Eq. 

The complete amplitude can now be written schematically as 

•Mtot = M-prod + M-tdecay + -M-tdecay ■ (6) 

The cross section, obtained from taking the absolute square of Mtot, then contains separate 
production and decay contributions, from |.M P rod| 2 and | M.tdecay\ 2 , \-M-tdecay\ 2 , respectively. 
It also contains cross terms representing the interferences, which in principle confound the 
separation but in practice are quite small. 

The interference terms are interesting in their own right, although not for top recon- 
struction. In particular, the interference between production- and decay-stage radiation can 
be sensitive to the top quark width T t 0, §], which is 1.42 GeV in the Standard Model at 
0(a s ) [[§. The interference between for example the two propagators shown in Eq. [| can 
be thought of as giving rise to two overlapping Breit-Wigner resonances. The peaks are 
separated roughly by the gluon energy, and each curve has width T t . Therefore when the 
gluon energy becomes comparable to the top width, the two Breit-Wigners overlap and there 
can be substantial interference. In constrast, if the gluon energy is much larger than T t , the 
overlap — and hence the interference — is negligible. Therefore the amount of interference 
serves as a measure of the top width. We will explore this more below. 

The integration over the final state phase space to obtain the cross section involves an 
integrand that contains multiple Breit-Wigner peaks from the top quark propagators as 

4 If the interference between processes is large this distinction is not even meaningful. 
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well as infrared singularities when the gluon energy becomes small. Even with cuts on E g , 
the rapid variation of the integrand can spoil the integration procedure. To eliminate this 
problem, we tailor the momentum generator to the production of a gluon in association with 
two massive particles (7*,Z* — > tig or t — > bWg). The multiple Breit-Wigner peaks are 
taken into account by using a multi-channel approach that integrates separately over the 
individual production and decay contributions; the Breit-Wigner behavior is smoothed out 
in the phase space generation. The interference terms, which have products of Breit-Wigners 
that peak in different places, much like in Eq. |], are are integrated using a combination of 
the three main channels. 

3 Numerical Results 

In this section we show results of the numerical calculation described above. We present 
the cross section for bbW + W~ g production in e + e~ collisions at a 500 GeV center-of-mass 
energy, with a few exceptions which are clearly identified. The calculation is entirely at 
the parton level, and we do not include initial state radiation, beam energy spread, or 
beamstrahlung. We use the following numerical values of parameters: m t = 175 GeV, 
nib = 5 GeV, Mw = 80 GeV, T t = 1.42 GeV, and a s = 0.1. Note that a s appears simply as 
an overall factor, because all of our events contain a gluon. 

Unless otherwise indicated, we use the following cuts. We require E g > 5 GeV to elim- 
inate the infrared singularity and because we intend for the gluon to be detectable. In 
addition we wish the gluon to be separable from the b and b quark; this is implemented via 
the requirement m,b g ,mi g > 10 GeV, which we shall identify below as "to& 9 cuts." (Separa- 
tion could also be achieved with a cut on the gluon's Et with respect to the b or 6; the choice 
makes little difference in the resulting distributions.) In order to make sure that we do not 
get contributions to our results from regions of the phase space where non-doubly-resonant 
diagrams might be important, we require 

160 GeV < m bW < 190 GeV or 160 GeV < m bWg < 190 GeV (7) 

and the same thing for the b. These conditions will be identified as "m t cuts" below. 

3.1 Characteristics of the gluon radiation 

We begin with the relative contributions of production- and decay-stage radiation to the 
total cross section. Figure |2] shows the fraction of the total cross section due to production 
stage emission in events with an extra gluon, as a function of the minimum energy of the 
gluon. This figure contains no cuts besides that for gluon energy and is simply meant to 
illustrate how radiation is apportioned in top production and decay for different center-of- 
mass energies; the solid line corresponds to cm. energy 1 TeV, and the dashed line is for 500 
GeV. Both curves fall off as the minimum gluon energy increases; this reflects the decrease 
in phase space for gluons radiated in the production stage. We see that the production 
fraction is always higher at 1 TeV collision energy than at 500 GeV. This too reflects phase 
space — for a given gluon energy there is more phase space available to produce gluons in 
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Figure 2: The fraction of gluon emissions radiated in the production stage, as a function of 
minimum gluon energy, for center-of-mass energy 1 TeV (solid line) and 500 GeV (dashed line), 
with no cuts. 

association with top pairs at the higher cm. energy. However both fractions remain below 
0.5; decay-stage radiation always dominates at these energies. 

Figure |] shows for a 500 GeV center-of-mass energy the effect on the production fraction 
of separation cuts between the gluon and b quarks. The dashed line shows the fraction with 
no cuts. The dotted line corresponds to requiring that the transverse energy of the gluon 
with respect to the b and b — which we denote Ex{g, b) — be greater than 3 GeV. The solid 
line corresponds to the cut m\> g > 10 GeV where the b can either be a quark or antiquark. 
The effect of both of these cuts is to eliminate gluons that are soft and/or close to one of 
the bottom quarks; since these contributions tend to come from decay-stage radiation, their 
effect is to increase the fraction of production-stage radiation. If the b were massless there 
would be a collinear singularity in the decay contribution; this does not happen in our case 
but the decay distribution still peaks when the 6-quark-gluon angle is small. The effects of 
both cuts become smaller with increasing gluon energy. 

Figure [| shows the total gluon energy spectrum for a collision energy of 500 GeV along 
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Figure 3: The fraction of gluon emissions radiated in the production stage, as a function of 
minimum gluon energy, for center-of-mass energy 500 GeV, with no cuts (dashed line), ET(g,b) > 
3 GeV (dotted line), and > 10 GeV (solid line). 

with its decomposition into production (dashed histogram) and decay (dotted histogram) 
contributions. The interferences between the two are negligible and are not shown; this will 
be true for all subsequent figures until we consider the interference explicitly. Included in 
this figure are the m t and ra& 9 cuts discussed above. As indicated in the previous figures, 
radiation from the top decays dominates. Otherwise the spectra are not vastly different; 
both exhibit the rise at low energies due to the infrared singularity characteristic of gluon 
emission, and both fall off at high energies as phase space runs out. 

3.2 Mass Reconstruction 

We now turn to the question of mass reconstruction. Figure |5| shows top invariant mass 
distributions with and without the extra gluon included; the first plot shows the distribution 
in rribw and the second shows nibwg- We have imposed rribg cuts and required E g > 5 GeV. 
In both cases there is a clear peak at the correct value of rrit. Note that the peak in the 
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Figure 4: The spectrum of radiated gluons as a function of gluon energy in GeV for center-of-mass 
energy 500 GeV, with m^g and mt cuts (see text). Dashed histogram: production-stage radiation. 
Dotted histogram: decay-stage radiation. Solid histogram: total. 

first plot contains the production contribution as expected, but the radiative decay part 
contributes as well. This is because even for decay-stage radiation, only one of the produced 
t quarks decays radiatively; the other still has p\ = p% w and therefore contributes to the 
rribw peak. The long tails in the two distributions are from misassignments of the gluons. In 
the left-hand plot, where the gluon is not included in the reconstruction, we see a low-side 
tail due to events where the gluon was radiated in the decay but was not included in the 
reconstruction. Similarly, in the right-hand plot we see a high-side tail due to events where 
the gluon was radiated in association with production, and was included when it should not 
have been. 

We wish to define a single distribution for the top mass that combines both types of 
events yet omits wrong combinations as much as possible. One possibility is to cut on the 
angle between the gluon and the b quark, whose distribution we show in Figure ||. This is 
motivated by the fact that gluons radiated from the b quarks are always part of the decay, 
and such gluons tend to be emitted close to the b quark direction. As we have mentioned, 
the mass of the b quark prevents a collinear singularity, but the gluon distribution still peaks 
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Figure 5: The top invariant mass spectrum without (left) and with (right) the gluon momentum 
included, for center-of-mass energy 500 GeV, with cuts and E g > 5GeV. 



close to the 6, as can be seen in the figure. Because we wish to define cuts that give a 
narrow top invariant mass distribution, the distribution in 9b g is decomposed into various 
invariant mass regions. (Here we refer to the b quark only, and not the b.) Using the variables 
fh t = m bW +, m t g = rribw+g, — m lw- an d fhf g = miw-g we define four types of events: 

type 1 : 172 GeV < fh tg , fht < 178 GeV (vertical hatching) 

type 2 : 172 GeV < fh t ,fht g < 178 GeV (horizontal hatching) 

type 3 : 172 GeV < fh t , fht < 178 GeV (cross hatching) 

type 4 : any other event (no hatching) 
Type 1 events are dominated by contributions from radiative t decays, and we can see that 
they do tend towards the b quark direction. Type 2 events (horizontal hatching) are in turn 
dominated by radiative t decays; gluons in this case tend to cluster near the b direction, 
and since the b and b tend to appear in opposite hemispheres, type 2 gluons are mostly 
found at large angles to the b. Events of type 3 (cross hatching) are mostly production-stage 
contributions; their distribution is more or less uniform, independent of the b quark direction. 
Finally, events of type 4 (no hatching) get contributions from both production and decay, 
with no compelling evidence for one over the other. 

Using this figure we can make the following conventions: 

if 6b g < 0.7 rad, assign gluon to t decay 

if 9i g < 0.7 rad, assign gluon to i decay 

if 6b g , 0b g > 1 rad, assign gluon to ti production. 
With these cuts on the proximity of the gluon to the b quark, we construct the top mass 
distribution presented in the dotted histogram in Figure [7]. 

Of course an important reason the cuts are so effective is that we work at the parton 
level. The experimentalists do not have that luxury, and, as one would expect, hadronization 
and detector effects are likely to cloud the picture. The solid histogram in Fig. |7| shows the 
mass distribution after including energy smearing; the solid curve is a Breit-Wigner fit. The 
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Figure 6: The distribution in the angle between the gluon and the b quark for center-of-mass 
energy 500 GeV, with E g > 5 GeV. The various contributions are as described in the text. 



spread in the measured momenta of the final state particles is parametrized by Gaussians 
with widths a = OAy/E for quarks and gluon, and a = 0.15\/~E for the Ws. We see that 
the central value does not shift, but the distribution becomes significantly wider. 

These results are meant to give an indication of the effects of hard gluon radiation on 
mass reconstruction and how they might be dealt with. Other variables to consider in 
choosing the cuts are m^, the transverse energy of the gluon with respect to the b or b, 
or some combination of energies and angles as defined in the various algorithms used in jet 
definitions for e + e~ colliders. At tree level and with partons only, the exact choice is not very 
important. We will revisit the question in more detail when we include virtual corrections 
in a full NLO calculation. 

3.3 Interference and Sensitivity to T t 

Finally, we return to the subject of interference. As mentioned above, the interference 
between the production- and decay-stage radiation can be substantial for gluon energies 
close to the total width of the top quark T t ; the interference is therefore sensitive to the 
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Figure 7: The top invariant mass spectrum with 6-gluon angle selection criteria (dotted his- 
togram), for center-of-mass energy 500 GeV, minimum gluon energy 5 GeV, and m^g cuts. The 
solid curve and histogram show the effects of energy smearing. 



value of T t . However, because this interference is in general small, we need to find regions 
of phase space where it is enhanced. This question was considered in || in the soft gluon 
approximation^, where it was found that the interference was enhanced when there was a 
large angular separation between the t quarks and their daughter 6's. 

Here we examine whether the result of |J, which considered a fixed final-state configura- 
tion, survives the exact calculation and phase space integration. Figure |S] shows that it does. 
There we plot the distribution in the angle between the emitted gluon and the top quark for 
gluon energies between 5 and 10 GeV and with cos 6 t b < 0.9 and m t cuts. The center-of-mass 
energy is 750 GeV. This cm. energy is chosen because for there to be significant interference 
between production and decay-stage radiation, both contributions must be sizable. At 500 
GeV, we see from Figs. and |] that the production contribution is suppressed compared to 
that from decay; as a result, the interference is very small. Increasing the energy increases 
the pro duct ion- stage contribution. We note that the distributions at 750 GeV and 1 TeV 
do not differ substantially. 



5 See also Q 10 
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Figure 8: The distribution in angle between the top quark and the gluon for gluon energies from 5 
to 10 GeV, cos t bi cos 8tb < ^-9, m t cuts, and 750 GeV collision energy. The upper solid histogram 
is the total and the other histograms represent the individual contributions: dotted: decay; dashed: 
production; dot-dashed: decay-decay interference; solid: production-decay interference. 

The histograms in Fig. |8| show the decomposition into the various contributions. The 
pro duct ion- stage radiation is shown as a dashed histogram; we see that it reaches its larges 
values at relatively small and large angles. Small angles correspond to the t direction, and 
large angles more or less to the i direction, since for the small gluon energies of interest 
here, the t and i are nearly back-to-back. The dotted histogram represents the decay-stage 
contribution; it dominates the cross section and peaks in the same region as the production 
contribution. This leads to substantial production-decay interference, shown as the negative 
solid histogram. This interference is destructive, so that it serves to suppress the total cross 
section, shown as the positive solid histogram. This effect would be enhanced if we lowered 
the gluon energies to values closer to T t , but jets from very low energy gluons are not likely 
to be observable, so we cut off the gluon energy at 5 GeV. Finally, interference between the 
emissions in the t and t decays are shown as a dot-dashed histogram, but as there is very 
little overlap between the two phase space regions even with these cuts, this contribution is 
negligible. 
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Figure 9: As in Fig. || with the addition of the cuts given in Eqs. P,[lO|. 

The cuts we have used are fairly generic; we can further enhance the interference terms 
with a judicious choice of additional cuts. If we examine their behavior in more detail 
in various regions of phase space, we find that the sign of the interference terms depends 
on the value of invariant mass of the top quark. Since we integrate over this mass, we get 
cancellations (a similar effect ensures cancellations of non-factorizable corrections in inclusive 
quantities). 

Consider the interference between radiation in the production stage and the top decay 
stage. The product of the two Breit-Wigner peaks is proportional to the factor 

ft = (Pwb ~ m t) * (PWbg ~ m t) + m t T t (8) 

This factor will multiply a quantity which, upon integration over angles, is negative. There- 
fore, for invariant mass values such that ft is positive, the interference terms are negative, 
while for negative ft, the interference terms are positive. We can impose cuts that take 
advantage of this: if we require the invariant masses to satisfy 

ft > , if e bg < 9- bg , (9) 

ft>o, ife bg >e- bg , (io) 
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Figure 10: The distribution in angle between the top quark and the gluon for gluon energies from 
5 to 10 GeV, cos Offr, cos Ojfo < 0.9, and 750 GeV collision energy. The histograms correspond to 
different values of the top width T t : dot-dashed: 0.1 GeV; solid: 1.42 GeV (SM); dashed: 5. GeV; 
dotted: 20 GeV. 

we obtain the distribution shown in Figure |9|. The interference effects are enhanced, though 
at the cost of a substantial decrease in cross section. 

Because the production-decay interference is destructive, increasing the top width would 
further suppress the total distribution. The height of the peaks, then, is sensitive to the 
value of T t . This is illustrated in Figure [It], which shows the cross section (without the cuts 
of Eqs. |9|JT0|) for different values of the top widthQ, ranging from 0.1 GeV to 20 GeV. The 
SM case (T t = 1.42 GeV) is shown as a solid line. It is interesting to note that in the con- 
text of perturbative gluon radiation, the SM top width is actually a small quantity. There 
are several other points to note. In principle, this sensitivity to T t gives us a method to 
measure the top quark's total width, independent of decay mode, above the top production 
threshold. Although in practice statistics would surely limit the possible precision of such a 

6 The histograms here are scaled so that they all would have the same normalization in the absence of 
interference effects. Without this rescaling, changing the width while keeping the tbW coupling fixed changes 
the total cross section, which behaves like l/T^ for small T t . 
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measurement, the total top width is not so easy to measure directly by any means. Further- 
more, the effects illustrated here arise from simple quantum-mechanical interference, and 
finding experimental evidence for interference between the radiation at the various stages is 
an interesting goal by itself. 

4 Conclusion 

In summary, we have presented results from an exact parton-level calculation of hard gluon 
radiation in off-shell top production and decay above threshold at lepton colliders, with the 
b quark mass and finite top width, as well as all spin correlations and interferences included. 
We have decomposed the cross section into the separate contributions from emissions associ- 
ated with top production, t and t decay, and their interferences. We have indicated some of 
the issues associated with this gluon radiation in top mass reconstruction and top width sen- 
sitivity in the gluon distribution. A detailed treatment of many of these issues, in particular 
with regard to mass reconstruction, will be presented in forthcoming work which combines 
virtual and real QCD corrections to this process into a complete NLO computation of top 
production and decay. 
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Appendix: Gauge Invariance 

Our results include only diagrams with two intermeediate top quarks. Because these top 
quarks are not assumed to be on shell, our result is not strictly gauge invariant. We can 
obtain a gauge invariant answer by subtracting some non-doubly resonant terms as follows. 
Consider the diagram where the gluon is radiated by the top quark (the first diagram in 
Fig. [I]). The contribution of this diagram to the production amplitude is (with a = a M 7 M ): 

M prod ~ qTZT" u (b) e w f f 6 + h (pt + k + m t )... v(b) . (11) 

zti Pt Pwb ~ m t 

By commuting e g to the right, this can be written 

M prod ~ wr- u{b) e w P J b + ™l (2e g ■ p Wb + e g k + (p 2 Wb - m 2 t )e g ) . . . v{b) . (12) 

* K Pt Pwb ~ rn t 

The term which breaks gauge invariance here is the one proportional to (piyt — rnf). This 
is a non-resonant term, regardless of the gluon being radiated in production or decay stage 
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(in other words, regardless of p\ Vh ~ m| or Pw bg ~ mf); therefore, in keeping with the 
approximation used, we can neglect it. 

A similar analysis works for the contribution of this diagram to the top decay amplitude, 
with py/b replaced by pwbg m this case, and we drop a term proportional to {p\ybg ~ m t)- 
Finally, the amplitudes corresponding to the diagram in which the gluon originates from the 
i can be computed in the same manner. The final result is gauge invariant, and differs from 
the exact result by non-doubly-resonant terms, as we have shown. 

We have implemented the above computation in the Monte Carlo program, and have 
checked numerically that the difference between the gauge invariant result and the exact 
result is very small (of order 0.01% of the total cross section, and order 1% with respect 
to the interference terms). This indicates that the other non-doubly resonant contributions 
(coming from diagrams with a single top or none) are also small; a more detailed study is in 
progress. 

Finally, we note that this method for restoring gauge invariance is not unique. We could, 
for example, have instead replaced the top mass in the top propagator numerator with the 
invariant masses: \jp\?\, in the production amplitude, and \Jpwb g m the decay amplitudes. 
The result is also gauge invariant, and also differs from the exact result by non-resonant 
terms. 
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